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We point out a new possibility for neutrinos where all neutrino flavors can be part Dirac and part 
Majorana. Our primary motivation comes from an attempt to use supersymmetric see-saw models 
to tie inflation, baryon asymmetry of the Universe and dark matter to the neutrino sector. The idea 
however could stand on its own, with or without supersymmetry. We present a realization of this 
possibility within an S3 family symmetry for neutrino masses, where we obtain tri-bi-maximal mixing 
for neutrinos to the leading order. The model predicts that for the case of inverted hierarchy, the 
lower limit on the neutrino mass measured in neutrinoless double beta decay experiments is about 
a factor of two larger than the usual Majorana case. 

I. INTRODUCTION 

Experiments over the past two decades have conclusively established that neutrinos have mass. The true nature of 
the neutrino mass however is unknown since available observations based on flavor oscillations do not tell us whether 
it is its own anti-particle (Majorana type) or not (Dirac type). Unlike the quarks and charged leptons, both these 
possibilities are allowed for the neutrinos since they are electrically neutral. Numerous neutrinoless double beta decay 
experiments are under way or in preparation to settle this question. 

An intermediate possibility that has been discussed in literature is known as the pseudo-Dirac case[l| where one 
includes a very tiny amount of the Majorana mass for each neutrino flavor which has Dirac type mass. The Majorana 
entry in this case must be very tiny (< 10~ 10 eV) in order to be consistent with current solar neutrino observations 0- 
In this note we point out a new class of possibilities where each neutrino flavor can have a large admixture of both Dirac 
and Majorana masses under certain circumstances. We point out the experimental implications of this possibility as 
well as its possible theoretical origin. 

While discussing the Dirac versus Majorana nature of neutrinos, it is usual to frame the discussion in terms of the 
neutrino flavor eigenstates that are emitted in beta decay or other weak interaction processes. When the neutrinos 
travel in free space, however, they do so as mass eigenstates, which are linear superpositions of the flavor eigenstates. 
This fact is responsible for neutrino oscillation phenomena. In this note we point out that the possibility of one 
of the neutrino mass eigenstates having a Dirac mass at the tree level with the others having Majorana type mass, 
appears consistent with all current observations. Since in this case, each flavor eigenstate is a large admixture of 
both Dirac and Majorana masses, we call this "schizophrenic neutrino" alternative. This is different from the usual 
pseudo-Dirac cases discussed in literature where the Majorana admixture is tiny compared to the Dirac mass whereas 
in our case Majorana admixture is as large or larger than the Dirac mass for each flavor. This possibility implies 
distinct predictions for neutrinoless double beta decay searches compared to the case where the neutrinos are pure 
Majorana type and could be used to test the schizophrenic hypothesis. 

On the theoretical side, the mass eigenstate having the Dirac mass must have a Dirac Yukawa coupling which is 
extremely tiny (~ 10~ 12 ) whereas the other masses could arise from high mass scale physics as in seesaw models [|[ 
with much larger Yukawa couplings. The Dirac type mass eigenstate would pair up with a right-handed (RH or sterile) 
neutrino (y s ) to form the Dirac mass. A priori, we do not know which of the three mass eigenstates has the Dirac 
nature. In this paper, we consider a specific model where we want to get tri-bi-maximal pattern [f| for the PMNS 
matrix. This suggests that the eigenstates be representations of an £3 symmetry group. The model then picks the 
"middle" eigenstate (the one that determines solar neutrino oscillations) as Dirac type since this is an S3 singlet 
whereas the other two are Majorana. In a generic model, any or in fact any two of the mass eigenstates could have 
Dirac type mass. 

While one would like to understand the small Dirac Yukawa coupling as a consequence of some high scale theory, 
it is comforting to know that it stable under radiative corrections due to chiral symmetry (or in this case under the 
symmetry v s —¥ —v s )- There may be other motivations for the existence of such tiny Yukawa couplings. One such 
motivation in supersymmetric versions of such models comes from attempts to use the RH sneutrino to drive infla- 
tion [J. In such a scenario, small Dirac coupling is essential to make the inflation predictions for density fluctuations 
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consistent with observations. However the hypothesis of schizophrenic neutrinos could be considered independently 
of this. As indicated earlier, a testable prediction of this model is that in the case of inverted hierarchy, the lower 
bound on mpp measured in neutrinoless double beta decay searches is roughly twice that of usual inverted hierarchy 
models in literature. This model will therefore be easier to rule out by the current generation of /3/3or/ experiments if 
long base line oscillation searches indicate inverted neutrino mass ordering. 

We hasten to point out that this kind of pattern for neutrino masses is not protected by a symmetry. As a result, 
when loop corrections are taken into account, tiny corrections of order g 2 m 2 /(32 v / 6tt 2 M2 / ) - 4 x 1CT 7 appear giving 
the Dirac eigenstate a pseudo-Dirac mass splitting of order 10~ 14 eV. These corrections have no impact on our 
prediction for Pffov decay. 



II. MOTIVATION FROM COSMOLOGY 



In this section, we review the cosmological motivation for small neutrino Dirac coupling in a supersymmetric seesaw 
model. We consider a supersymmetric extension of MSSM based on the gauge group SU(2) L x U(l)i 3R x U(1)b-l 
which requires that there be three RH neutrinos N (= v c in SUSY language) to cancel anomalies (with the eventual 
possibility of SO(10) grand unification). A combination of superpartners </> = (N + H u + L/^/3) in this theory is a 
Z?-flat directions under the whole gauge symmetry , and is also i^-flat when neutrino Yukawa couplings h are turned 
off. As shown in this fiat direction can act as the inflaton. The neutrino Yukawa couplings in combination with 
the soft mass and A-term for <p leads to a potential for the radial component of <fi (denoted as if) of the form 0] 

v{ v ) = + T y - (i) 

where m 2 = (m 2 ^ + rn 2 Hu + m|)/3 can lead to inflection point inflation 4] and the amplitude of observationally 

relevant density perturbations (as measured by COBE and WMAP) matches the observed value Sh ~ 1-9 x 1CP 5 for 
weak scale masses ~ 0(100) GeV, provided that h ~ 10~ 12 (for details, see 0). The neutrino mass is intimately 
connected to h. For instance, if neutrinos are Dirac fermions, we have m„ = h(H u ), where (H u ) = (174 GeV) sin/3 
and tan/3 is the ratio of vacuum expectation value (VEV) of Higgs fields of the minimal supersymmetric standard 
model (MSSM). Then h ~ 10~ 12 would give rise to m u ~ O(0. 1) eV, which is precisely in the range of interest for 
neutrino oscillations. We could take this as a hidden message from cosmology that at least one of the neutrinos can 
be dominantly of Dirac nature, and study its implications for neutrino masses and mixings. 

It is important to emphasize that the inflation model constrains the coupling of only one of the RH neutrinos whose 
superpartner is responsible for inflation. That RH neutrino could be the Dirac partner of one of the light neutrino 
combinations making it a Dirac neutrino. The other two RH neutrinos have unconstrained Yukawa couplings that 
take natural values (~ 10 -5 — 0.1), and hence their mass must be heavy. Note however that the heavy RH neutrinos 
must not mix with the RH neutrino whose superpratner is part of the inflaton so as not to spoil the picture of inflation 
mentioned above. The simplest possibility for neutrino masses in this case would therefore appear to be that one 
linear combination of the flavor eigenstates is a Dirac fermion whereas the other two will be Majorana and get their 
mass via the see-saw mechanism. Below we suggest this as new picture for neutrino masses. 



III. AN S 3 MODEL FOR SCHIZOPHRENIC NEUTRINOS 



One of the challenges in neutrino mass physics is to understand the observed near tri-bi-maximal mixing pattern 
among different flavors. Discrete symmetries have been discussed extensively as a way to address this issue @ and 
the group S3 is one of the symmetries that appears promising in this context and we use it in our discussion in this 
paper. The basic assumptions of our neutrino model can therefore be summarized as follows: 

• The extended gauge group responsible for neutrino masses consists of a local B — L symmetry 0? ], which 
requires the existence of three RH neutrinos for anomaly cancellation. 

• One of the RH neutrinos (whose superpartner is the inflaton field) couples to a linear combination of all neutrino 
flavors with a Yukawa coupling of order 10 -12 so that it gets a Dirac mass without any need for see-saw, whereas 
the remaining orthogonal combinations get their masses from the conventional see-saw mechanism. 



• The three standard model lepton doublets transform into one another under a flavor S3 discrete symmetry. 
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The first assumption is quite well motivated and has been widely discussed in literature. It also naturally incor- 
porates N along with H u and L into a single D-flat direction that can drive inflection point inflation. The second 
assumptions is motivated by the cosmological scenario discussed above. 

As already mentioned, our neutrino model is based on the idea that only one of the neutrino flavor combinations 
corresponding to a mass eigenstate has a small Yukawa coupling to one RH neutrino whereas the other two combina- 
tions get their masses from the seesaw mechanism. If we take the tri-bi-maximal matrix as the leading order PMNS 
matrix, then one might start thinking of a discrete symmetry group which has one singlet and one doublet as part of its 
irreducible representations and the singlet one being the Dirac neutrino whereas the doublet combinations becoming 
Majorana. One such example used in literature is the S3 group [Ij which proves convenient for our discussion. 

We assume the S3 to permute the three families of leptons (£ e , L^, Zr T ) among themselves. Of course, it is well 
known that this is a reducible representation of S3 group but the following linear combinations of these fields transform 
as singlet and two dimensional representations of 5*3: 

L 2 = -j=(L e + + L T ) (Singlet) 

Li = -=(2L e - Lp- L T ) (Doublet) 
V6 
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L 3 = —(L^-Lt) (Doublet). (2) 



We assume that muon type RH neutrino is the S3 singlet whereas (N e ,N T ) form a doublet. The masses of these 
doublet RH neutrinos can be chosen different by appropriate choice of symmetry breaking (see below) . The effective 
lepton Yukawa coupling after integrating out N e and N T can then be written as 



C v = hL 2 H u N ll + -^-(L^) 2 + ^-(L 3 H u ) 2 

M Nc M Nt 



+ h.c. (3) 

After the electroweak symmetry breaking, this gives rise to one Dirac neutrino corresponding to the mass eigenstate 
V2 and two Majorana eigenstates ^1,^3 and clearly leads to tri-bi-maximal form for the PMNS matrix provided the 
charged lepton mass matrix is diagonal. 

The effective Lagrangian in [3] could for instance arise in an SU(2)l x U(1)i 3R x U(1)b-l theory supplemented by 
a global discrete symmetry S3 x D, (where D is a product of extra Z n symmetries) if we take an S3 doublet Higgs 
fields (Af , Aj) and a singlet field Ag (all with B — L charge —2 and I3R charge +1) with VEVs (Ai) = and others 
with non-zero VEV. This will generate different Majorana masses and M^ T for the S3 doublet RH neutrinos. 

In this model, inflation occurs along the flat direction corresponding to the first superpotential term in Eq. [3] The 
coupling between and (N e , N T ) can be forbidden by e.g., a Z% symmetry contained in D under which —iN^ 
and a gauge singlet field X with X — > e I7r / 4 A with all other fields invariant. The Dirac coupling of can be obtained 
from a higher dimensional coupling (XL2H U N I1 X 2 /Mj,,), where (X) ~ 10 12 GeV or so. At the inflection point VEV 
(~ 10 12 GeV [I)]), this interaction then yields the effective Dirac coupling of in Eq. |3J An additional RH neutrino 
mixing term (NnN eT A° 2 A 2 /Af|, i ) is allowed by the Zg symmetry, but has negligible contribution to masses and 

mixings and can be ignored. We need to add the fields (A^, A 2 ) and a singlet field Aq to preserve supersymmetry 
below the B — L scale as well as to cancel anomalies. 

Turning to the charged lepton mass matrix, neutrino oscillation data require that it be nearly diagonal. We employ 
the technique used in the second reference in Q. We add three gauge singlet superfields (a e , a^, a T ) and three extra 
Z n symmetries, i.e. Z n>e x Z n4L x Z n>T , with RH lepton fields e c ,^ c ,r c transforming as and singlet fields 

transforming as T . Both sets of three fields also transform under S3 exactly like the lepton doublet fields. We 
can write down the corresponding Yukawa superpotential as 



Wiy = JT^Hd {L e a e e c + L^o^f + L t <t t t c ). (4) 



There can be another term where the L e ,L^,L T are permuted among themselves. This will contribute to the off- 
diagonal elements of the charged lepton mass matrix after symmetry breaking. We set this coupling to zero. Now 
by adjusting the VEVs of the singlet fields, we can get diagonal mass matrix for the charged leptons. On the other 
hand, if the small contributions to the mass matrix coming from the permuted terms are kept, there will corrections 
to the tri-bi-maximal form e.g. it will lead to non-zero #13. 
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FIG. 1: We plot Im^l as a function of the lightest neutrino mass for the case of inverted hierarchy. The red (dark) band shows 
the prediction from the two Majorana and one Dirac neutrino scenario and the gray shaded region shows the conventional three 
Majorana neutrino case. 



IV. IMPLICATIONS FOR NEUTRINOLESS DOUBLE BETA DECAY 



This neutrino mass model has an interesting implication for neutrinoless double beta decay. Recall that in the con- 
ventional all Majorana neutrino case, the light neutrino contribution to /3/3o„ decay takes the form mpp = U^rrii, 
(i = 1, 2, 3), where U e i are entries of the PMNS matrix. In the inverted hierarchy scenario, this leads to the following 
lower bound for the conventional three Majorana neutrino case Q 

\m pp \ ~ 1(008% + e M sin 2 )m atm | > ^ « 17 meV 

(Conventional). (5) 

In our model, however, the second neutrino mass eigenstate is a Dirac type state and therefore has no contribution 
to PPov decay. This leads to the following lower bound for inverted case: 

\mpp\ ~ cos 2 6» m atm > 2m * tm « 34 me V 

(Dual), (6) 

which is roughly twice the value of the conventional case ([5]). This makes it easier to rule out our model in the current 
generation of neutrinoless double beta decay searches, provided we have independent evidence, e.g. from long base 
line neutrino experiments for inverted mass hierarchy. 

In the normal hierarchy scenario, the corresponding formula becomes mpp ~ (miCOS 2 #o+e IQ sin 2 0i3m atm ), which is 
different from the conventional three Majorana case. The precise value in this case depends on the unknown Majorana 
mass of v\ as well as the value of #13. 

In Fig. [TJ we plot Im^l as a function of the lightest neutrino mass min(mj) = m m i n (which sets the absolute 
neutrino mass scale in the case of inverted hierarchy). The red (dark) band shows the prediction of our scenario 
and the gray shaded region shows the usual three Majorana neutrino scenario in the case of inverted hierarchy. The 
masses and mixing angles used for the figure are as follows [10]: Atoq=(7.59 ± 0.20j^Qgg) x 10~ 5 eV 2 , Am 2 tm = 

(2.46±0.12±0.37) x 10- 3 eV 2 , 9 13 < 12.5°, Q = 34.4° ± l^; 2 and 9 atm = 42.8° t^'l-W ■ We can see that in the case 
of inverted hierarchy (corresponding to m m i n < 0.05 eV) the lower limit on |mg^| measured in neutrinoless double 
beta decay experiments is about a factor of two larger than the conventional case. 



V. COMMENTS 



We now make some comments on the model described here. 

(i) Since the Dirac nature of the second neutrino mass eigenstate is not protected by any symmetry, radiative 
corrections will induce Majorana component to its mass. The self-energy corrections to Vi masses due to W + £~ 
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intermediate states will lead to kinetic mixings between the different mass eigenstates that depends on the charged 
lepton masses: ~ Cl2k UikU£ g 2 m 2 k 32^/6TT 2 M^) . This mixing is of order 4 x 10~ 7 . When the kinetic energy term 
in the Lagrangian is diagonalized, this leads to mixing terms (for the normal hierarchy case), e.g. msv^V2 + where 
ma ~ TO^ 3 e23 + This introduces a Majorana mass term hm-iv-iv-i with Sm,2 ~ 10 -14 eV. It, being very small, keeps 
the Dirac nature of v-i to very high precision. This is also consistent with a bound ~ 10 -9 eV on this parameter from 
solar neutrino observations |2|. The same result holds for the inverted hiercrhy case with v\ and v% interchanged. In 
the SUSY version, quantum corrections that mix the slepton states introduce a Majorana component for the Dirac 
neutrino. At one loop this effect results in Smfj ~ [(Y^Yi)^ /Ww 2 ] TOq ln(M 2 /Af§), which is of the same order as 
that mentioned before. 

(ii) We wish to emphasize that our scenario is different from the usual pseudo-Dirac scenario [l| discussed in the 
literature. Our light neutrino mass matrix is a 4 x 4 matrix such that one of its eigenstates forms a Dirac pair with 
the sterile neutrino and other two eigenstates are Majorana. The Dirac eigenstate gives rise to a pseudo-Dirac pair 
only at the one-loop level. 

(iii) The masses of the two heavy RH neutrinos depend on the scale at which B — L is broken, and can be as low 
as 0(1) TeV. The decay of heavy Majorana neutrinos, and their SUSY partners, can generate baryon asymmetry of 
the Universe. If Mtv c ,Mjv t are of order TeV, resonant leptogenesis will be a relevant solution. However, in the S3 
symmetric model, this docs not work since it will require the first and the third neutrino masses be almost equal. The 
oscillation data will be hard to fit with this pattern. However, soft leptogenesis [ll| can work well in the model for a 
wide range of Majorana masses. 

(iv) In this model, either the MSSM neutralino or the superpartner of the RH component of the Dirac neutrino 
can play the role of dark matter. The latter is naturally the lightest of the RH sneutrinos since its mass receives 
contribution from SUSY breaking alone. If the B — L is broken around TeV, it can obtain the correct relic density 
via thermal freeze out [4|. This also makes the corresponding Z' accessible at the LHC. On the other hand, for a high 
scale B — L the usual MSSM neutralino is a good dark matter candidate. The role of B-L in this case is to provide 
the R-parity symmetry naturally. 

(v) The impact of the RH neutrino, which is responsible for the Dirac mass, on Big Bang Nucleosynthesis also 
depends on the scale at which B — L is broken. For example, for Mz> ~ 10 TeV, the RH neutrinos decouple at 
T D ~ 1 GeV, while for Mz> ~ 1 TeV we have T D ~ 100 MeV. In the latter case this amounts to ~ 4, whereas 
N° s ~ 3.1 in the former case. For a high scale B — L the RH neutrinos decouple much earlier, and hence N° s w 3. 

In summary, motivated by cosmology, we have pointed out a new picture for neutrino masses with the novel 
feature that one of the mass eigenstates is a Dirac fermion (at the tree-level) whereas the other two are Majorana 
type. We presented an S3 realization of this idea that leads to tri-bi-maximal mixing for leptons in the leading order. 
This model can be ruled out by the current generation of neutrinoless double beta decay searches if inverted mass 
hierarchy is indicated by long base line neutrino oscillation experiments and neutrinoless double beta decay searches 
give 17710/3 1 ^ 34 meV. 
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